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I. INTRODUCTION 



Random matrices are known to find applications in many physical systems [H, Q • I n the present work, we focus on 
rotationally invariant ensembles of N x N matrices H (symmetric, hermitian or quaternion self-dual), for which the 
joint probability density (jpd) of the N real eigenvalues {Ti} can be generically written as: 



P(Tr,..., T N ) = -—— J] \Tj T k f J] u (Ti) (1) 
^^ iy ) j<k i= i 

where (3 is the Dyson index of the ensemble (f3 = 1, 2,4 respectively), ui(x) a certain weight function and Z w (/J, iV) 

is the normalization constant. The classical ensembles of random matrix theory correspond to the following weight 
functions: 

2 /o 

cu(x) = e~ x * — oo < x < oo Gaussian Ensemble (GXE) (2) 

u)(x) = e~ x x a ~ 1 x > Laguerre Ensemble (LXE) (3) 

Lu{x)=x a - 1 (l~xy- 1 0<x<l Jacobi Ensemble (JXE) (4) 

where X = {0, U,S} stands for Orthogonal, Unitary and Symplcctic ((3 — 1,2,4 respectively). The normalization 
constant Z U ([3,N) can be computed via the celebrated Selberg's integral: 

„i „1 n 

S n (a,b,c):= •••/ dt 1 ---dt n \[t a l - 1 (l-t i ) b - 1 J] l*i-*j| 2c 

" _1 T(a+ jc)T(b + jc)T(l + (j + l)c) 
T(aH 

and its generalizations. In particular, we have 



=o T(a + b+(n + j-l)c)r(c + l) {) 



^ T f 1 + ij) 

Zueoifi, N) = (2nf/* J] > J (6) 



j 

N-l 



L ir(i + f) 



7 IR n\ TT r(a + j/3/2)r((j + l)/3/2) 
Z^j(P,N)=S N (a,^0/2) (8) 
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In many physical applications, one is interested in so-called linear statistics on the N eigenvalues, i.e. random 
variables of the form: 

N 

A = J2f(T>) (9) 

8=1 

where the function f(x) may well be highly non- linear (see 0, 0, H, [|| and references therein for physical applications). 
In particular, general methods are available [1,0] to compute in principle the mean and variance of any linear statistics 
from a generic invariant ensemble, at least in the large N limit. 

Conversely, much less is known for nonlinear statistics, i.e. functions involving products of different eigenvalues 
(see however [8, 9]). One may for example consider the following random variable: 

N 

T* = perm(¥) := £ \{^(i){Ti) (10) 

where perm stands for the permanent of the N x N matrix \I/ = ('0i(7 1 j))i<j.j<A r , the sum runs over the permutations 
7r of the first N integers {&n is the symmetric group), and {ipi(x)} is a set of iV given functions. Clearly, the general 
definition above: 

1. is invariant under permutations of the Tf, 

2. incorporates as special cases e.g. powers of the determinant of H, ((det H) K ) [l(| |Tl[ [eq. (|10P when tpi(x) — x K Vi] 
as well as traces of higher powers of H, (TrH K ) [cq. (JTDJ) when ipi(x) = x K for i = 1 and 1 otherwise]. 

The aim of this paper is to study the statistics of permanents on classical random matrix ensembles. Our 
motivation comes from the problem of quantum transport in open chaotic cavities supporting N\ and N2 electronic 
channels in the two attached leads. A detailed account of the problem and its link with the Jacobi ensemble of random 
matrices is provided in Appendix [AJ Our more general approach allows to extend results from [1 21 ] in a clear and 
computationally efficient way. 

We base our analysis on the theories of hyperdeterminants and symmetric functions. The former is detailed in 
Appendix [Bl and involves multidimensional generalizations of the conventional determinant: it provides a very general 
(although not always efficient) way to write averages of permanents as sums of determinants for (3 — 2. The latter 
is detailed in Appendix [C] and will be used to produce less general, but quite powerful formulae for a few physically 
interesting cases. The use of hyperdeterminants and symmetric functions in random matrix contexts is not new (see 
e.g. [H, 0, [H, [H, [TtJ and references therein). Here we apply similar methods to a different problem. Note that 
the links between Selberg integrals and hyperdeterminants have been already investigated by one of the author with 
Jean- Yves Thibon QJ, Gj| . 

The plan of this paper is as follows. In section [TT] we provide a general hyperdeterminant formula for the average of 
permanents valid for (3 = 2, any St and any benign weight lu(x). While being very general, the practical implementation 
becomes rapidly unwieldly due to an exponential growth of the number of terms with N. In section ITOl we resolve this 
efficiency issue adopting symmetric functions expansions. The resulting formulae do not increase in complexity when 
N grows, thus making the numerical implementation extremely efficient, though at the price of a loss in generality. 
These improved formulae are valid for (3 = 1,2,4 for the Jacobi weight with 7 = 1 and we restrict ourselves to the 
most interesting case ipi( x ) — xXi ■ Finally in section HVl we offer some concluding remarks. In the appendices, we 
give a detailed account of the problem of quantum transport in chaotic cavities which constitutes our motivation (|A|) 
and some remarks abour hyperdeterminants JB]) and symmetric functions |[C[). In (|Dj) . guided by the numerics, we 
perform an asymptotic analysis for large N and put forward a factorization conjecture that will be studied in more 
details in a forthcoming paper [20j . 



II. HYPERDETERMINANT FORMULA FOR STATISTICS OF PERMANENTS FOR (3 = 2 

We are now interested in computing the following average: 

1 r N 

<perm(*)} = / J] dT, uft) ]J \T, - T/perm^fT^K,,^ (11) 

where uj[x) is one of the classical random matrix weights and the integrals run over the appropriate support. Hereafter 
it is assumed that both the measure u>(x) and the functions il>i{x) are benign, i.e. they ensure existence and convergence 
of the integrals involved. 



3 



Example 1 Consider tpi(x) = x Xi , where A = [Ai, A2, •• • , Xn] is a decreasing partition (Ai > A2 > . . . > Xn). Then: 

(perm(*))=A! (T^---T^) (12) 

Definition (fTTj) is very general: it requires to specify the Dyson index j3, the measure u>(x) and the set of functions 
{ipi(x)}. In this section, we will focus mainly on the unitary case ((3 = 2), all the remaining 'degrees of freedom' 
being left untouched. In principle, the same reasoning could be applied to any even /?, but the resulting formulae are 
too complicated for any practical use. 

The main technical tools are the following: 

• The expansion of a hyperdeterminant as a sum of conventional determinants (see eq. (|B4[) in Appendix [Bj . 

• A generalization of Heine's theorem for determinants 1 . Given £ sets of A functions ft (x) (s = 1,...,£ and 
i = 1, . . . , A) and a benign integration measure n(x) the following holds: 

N k I 

JJe^(xi) Y[peTm(f^\xj))t<i tj < N det(tf s) (xj))x<i,j< N = 

= 1 s=l s=fc+l 

= A! Det {fc+1 ,... /} ( f dfx(x) J] /« (x) j (13) 



l<ii,...,V<iV 



• The Cauchy's double alternant evaluation [23|. For any pair of vectors X = {Xx, . . . , X n } and Y = {Yx, . . . , Y n } 
the following holds 



1 \ IllKKjKn&i-XjHYi-Yj) 



det --fy = (x+Y) (W) 

\A l + *jj x<i,j<n \.\.l<i,j<n\ A i + y ii 

Let us introduce for (3 = 2 a special case of hyperdeterminant Det+ (see the general definition and properties in 
Appendix [B| of the multi-indices tensor M defined by 

1 N 

T>et + (Mij, k )x<i,j,k<N = jfi e ( <J 20- 3 ) J|M CTlWi0 . 2(i ) )(T 3 (i) . (15) 

0-1,0-2, o- 3 es^ »=i 

where the sum runs over permutations ax, o~2, 03 of the first N integers and e is the product of their signatures. 
In the following, we will make use of two of the properties just stated above, namely: 

1. The expansion in terms of conventional determinants (see again Appendix |B|) 

Det + (M <iii fe)i< i j i fe<iv = Y det ( Af o(i),i,i)- (16) 

ctE&n 

2. The generalization of the Heine theorem for hyperdeterminants (fT3"|) : 



JV 

perm(fi(xj))x<i,j<N det(gi(xj))x<i,j<N det(hi(xj))x<i,j<N u(xj)dxj 



A! Det+ / dx u(x)f % {x) 9j {x)h k {x) . (17) 



l<i,j,k<N 

From Heine's theorem, the evaluation of the average pT|) for /3 = 2 is quite straightforward. Noting that: 

n(T,-T fc ) 2 =det(77'- 1 )^. ij < Ar (18) 

j<k 



1 The original Heine's theorem can be found in 1211 . Conversely, eq. I)13|l does not appear explicitly in literature but a version for totaly 
alternated hyperdeterminants can be found in |22l and its proof is straightforward. 
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one simply has: 



(perm(*)> - 7 ^ r^ Det + f / dx ^{x)^)^^ 2 ] (19) 

Z„{Z,1\) \J / Ki,j,k<N 



l<i,j,k<N 

From the first property stated above, this can be expanded as a sum over the symmetric group & n 



(perm(*)) = Y~JYW) ^ ^ (/ ^ 



(20) 



Eq. (|20p is the main result of this section. It expresses in full generality the average of any permanent for any 
invariant ensemble with /3 = 2 as a sum of N\ determinants. Clearly, due to the exponential growth with N of its 
complexity, formula (|20|) is only practical when N < 10. We will offer in section 11111 a less general, but more powerful 
way to compute the sought average. 

Example 2 Suppose we take the Jacobi weight with 7 = 1, lo(x) — x a ~ l (relevant for the quantum transport problem) 
and the nonlinear statistics il>i(x) = x Xi , where A = [Aj., .. . , Xn] with Ai > ■ • • > Ajv- Then formula (|20p reads: 

/ T Ax T x N] _ (perm(tt)) _ ^e6 ff det (jp^ a ' a+A " i)+ ' +, " 3 ) 1 ^ gjV _ ^,ee N det ( a+K J +i+j - 2 ) ^. ^ 

(21) 

Further simplifications are achieved employing the Cauchy 's double alternant identity, which finally yields: 

/mXi t X n k _ UiKji 1 - j) sr-y I[\<i<j<N^(i) + 1 ~ ~ j) . , 

Z^J(2, N) ^ n£=i(A«r(i) + * + J + « - 2) 



(|22p extends in a compact form the results from Jla] to arbitrary values of Aj. In the special case A = (1, . . . , 1), 
we note that (|22| perfectly matches the value for the average of the determinant for a Jacobi ensemble with 7=1, as 
computed from the Selberg integral (J5|) : 

(det(H)) Q = (Ti • • • Tn)u = S f ( +1 :]f = n (23) 

b N {a,l,l) a + N + j 

thanks to the following (easy to prove) identity valid Va, n: 

n] Ui^zJl = rf r (" + J + i) r 0' + i) r 0' + 2) f24 x 

nk=i(*+i+«-i) ;i r(a+i+n+j) 1 ; 



III. MORE EFFICIENT SYMMETRIC FUNCTION EXPANSIONS FOR THE JACOBI WEIGHT 

In this section, we are able to provide more efficient and user-friendly formulae for a special nonlinear statistics, 
namely the quantity: 

1 r N 

(T^ —T*«) a = — — / dT x ■ ■ ■ dT N T x Al ■■■T^Yl \Tj -T k f]J T^ 1 (25) 

where the average is taken with respect to the Jacobi weight @ with 7 = 1 and /3 = 1, 2, 4. 

Such object is of interest for the statistics of moments of experimental observables in the problem of quantum 
transport in ballistic chaotic cavities. A detailed overview of the problem is provided in Appendix [XJ It is likely that 
the method we present here, based on symmetric function expansions, may be applied with slight modifications to a 
number of other measures and observables. 

The main tools are the following: 
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• The following identity (hereafter we use the notation of [24]): 

perm(Tf 3 )i<ij<N = \ ] m\ (26) 
where m\ is the monomial symmetric function (24| 

m x := mx(Tx, . . . , T N ) = ^ r/ 1 • • • (27) 

i 

summed over all distinct permutations J of A, and X' = uqI . . .n^- ■ ■ ■ if rij denotes the number of occurences of i in 
A (for example, [5, 5, 5, 3, 3, 2, 1, 0, 0, 0] ! = 3!0!2!l!l!3!). 

• The well-known link between Selberg-type integrals and Jack polynomials given by the Kadell formula [25| : 

„ N 

jj? := / j[°HTi,...,T N )]l\T i -T j \ 2c mi-T i ) b -« N - 1 ^^^ 

_ 1 r(d + i)r(b + c{i + i))r(A. t + a + c (i - p) 

~ Ja r(c + i)r(A i + a + 6 + c (i- ? )) • (28) 



1=1 



(I) 

where the value of J\ " (1, . . . , 1) is known to be 

jyp.-.D-^n r^-.+i)) (29 » 

if A^ > ^(A) (where £(X) denotes the length of the partition A and |A| the sum of its nonzero parts) and otherwise 
(24J. The Jack polynomials for c = 1 are proportional to Schur functions and for c = 1/2 are known as zonal 
polynomials (see Appendix [Ul for details). 

We first illustrate in some detail the method for (3 = 2 (c = 1) based on Schur function expansions, where we 
provide also a detailed asymptotic analysis for N — > oo, and then the cases j3 = 1,4 together on the same footing. 



A. The Unitary case (/3 = 2) 
Consider the expansion of m\ in the Schur basis [24j 

m x = Y,KW- (3°) 

The coefficients are obtained inverting the Kostka matrix [24| . an operation that can easily be performed by 
symbolic computation routines. 

Replacing the permanent by its expansion in the Schur basis one finds (for a = 1): 

\! ~ f N 

(Tt...T^) a=1 = -— J2 K x s fl (T ll ... 1 T N )det(Tr 1 )i <l , J<N l[dT 3 (31) 

A! Z u =j{2, A) J [0>1] n -' J - 

The integral on the right-hand side is readily recognized as a special case of Kadell's integral (|28[) for c = 1 , leading 
immediately (after simplifications) to the final result (eq. ([3~4")> and (J35J). However, it is more instructive to proceed 
directly from (|3~lj) and notice that each Schur function is itself the quotient of two determinants [2~4| : 



d et(r I ^ +JV - J ) 1 < M -< w 

de^T/" 1 )^^^ 

It follows that the integral in (|3Tj) becomes: 



^m, . . . , t n) = uc ^;^._ (32) 



(T^ ...T A «) Q=1 = -— A ' ^ E^a / <!<'•;/;' ' V W v-l.'-i/r x\['IT, (33) 

A I Z uj =j(2,1\) ^ J[o,i]H ^i 
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and using the Heine theorem, each multiple integral in the sum can be again converted into a determinant. This 
procedure can be easily implemented for weight functions different from Jacobi, and will lead to general expressions 
for averages like (T^ 1 . . - T N N ) as computationally efficient sums of determinants. In the present case, following one 
strategy or another and exploiting the Cauchy's identity (TT4")) we obtain: 



The parameter a can be introduced easily: 



n 



! 1 1«<J 



(i ~ 0) n i <, ifJ-i - IH +3-i) 



(2, AO 



(T^...T^) a =\ 



(i-3) 



Z usJ (2,N) 



(34) 



(35) 



The formula (|35p is the main result of this section. It provides a very efficient algorithm (see examples below) , since 
the size of the sum does not depend on the size N of the alphabet (compare it with formula (j22|) ). Knowing the Schur 
expansion of the monomial, the computation is immediate. Let us illustrate this method on the following example. 

Example 3 Let X = [4, 3, 2] . One has 

m [4,3,2] = s [4,3,2] - s [4,3,l,l] ~ S [4,2,2,l] + 2 S[4,2, 1,1,1] 

— 2 S[4, 1,1, 1,1,1] — 2S[ 3 ;3 3 ] + 5(3,3,2,1] — 2S[ 3 ,2,i,i,i,i] 
+4*[3,1,1,1,1,1,1] + 2 S[2, 2, 1,1, 1,1,1] 
-6S[ 2 , 1,1, 1,1, 1,1,1] + 6S[i, 1,1, 1,1, 1,1, 1,1] 



Hence, 



(TfT^T^a — /[4,3,2] — /[4,3,1,1] ~~ /[4,2,2,1] + 2/(4,2,1,1,1] 

— 2 /[4, 1,1, 1,1,1] — 2 /[3,3,3] + /[3,3,2,1] ~ 2 /[3,2, 1,1,1,1] 
+4 /[3, 1,1, 1,1, 1,1] + 2 /[2, 2, 1,1, 1,1,1] 

— 6 /[2, 1,1, i,i, 1,1,1] + 6 /[i, 1,1, 1,1, 1,1, 1,1] 



where 



/u = A ! 



n i<7 (* - 3) il< (Mi - ftj +3-*) 



Z u =j(2, N) Ui,M +N-t+j + a-l) 
if the length of fi is less or equal to the number N of variables Tj and otherwise. 
• For N = 3, one finds after simplifications 



(T?TiTi) a 



(28 + 10 a + a 2 ) (2 + a) 2 (3 + a) a (1 + af 
(6 + a) 2 (4 + a) (7 + a) (5 + af (8 + a) 



• For N = 8 



with 



• For N 



with 



P(a) (5 + a) (6 + a) (7 + a) 



(11 + a) (9 + a) (12 + a) (10 + a) (13 + a) (14 + aT (15 + a) 3 (16 + a) (17 + a) (18 + a) 



P(a) = 422781389568 a + 166843216800 a 2 + 40063436856 a 3 + 6512032020 a 4 + 3924093 a 7 + 63580545 a 6 
+753772094 a 5 + a 11 + 105 a 10 + 174690 a 8 + 5388 a 9 + 493650339840 

20 we have (after a computation of few seconds on a standard personal computer) : 

r4rr3rr2 , _ Q(a)(17 + a)(18 + a)(19 + a) 



{T?T$T£) C 



R(a) 



Q(a) = 17973994269257913600a + 3006012942356996160a 2 + 311313563528661024a 3 
+1371214528697a 7 + a 13 + 360a 12 + 45753367235370a 6 + 1164250996862956a 5 + 
62879a 11 + 6977370a 10 + 31407665820a 8 + 544626843a 9 + 50279359153701888000 



22279414065220920a 4 



• For N 



R(a) = (31 + a) (32 + a) (33 + a) (34 + a) (35 + a) (36 + a) (37 + a) (38 + a) 2 (39 + a) 3 (40 + a) 2 (41 + a) (43 + a). 
3 ... 50, see Fig. LH 



Asymptotic analysis of the unitary case for N — » oc. This analysis is reported in Appendix ID] 
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FIG. 1: Values of {TtT'iTi) a as a function of a € [0,50] for /3 = 2. From (jD6j and fDT5| for p < 1, one has 
lim J v_ o<r 1 4 r 2 3 r 3 2 ) a = 525/16384 ~ 0.032 in good agreement with the plot. 



B. The Orthogonal (f3 = 1) and Symplectic case (/3 = 4) 

In complete analogy to the case = 2, the algorithm to compute (very quickly) (T^ 1 . . . Tjy ) a from (|25|) consists 
of three steps: 

1. Replace the permanent (l/Af!)perm(Tj A3 )i<ij<jv = 7\ Al ■ ■ - T^ N with a monomial symmetric function using (|26)l ; 



2. Expand the monomial function m\ in the Jack basis for the parameter c = 2 (/3 = 4) or c = 1/2 (/3 = 1) [24j; 

( 2 ) ! ( 2 ) 

3. Replace each occurence of by jjjZu=j{(3, N) -1 !^ 
Let us provide a couple of examples for such procedure: 

Example 4 Consider the average (TfT^T^a for (3 = 4. The expansion of the monomial function 771(4,3,2] in the Jack 
basis is (for an alphabet of size N = 3) 

m [4 ,3,2] = -j^J} 3 %] + 2-^ J WY ( 36 ) 



After substitutions, one obtains: 



/■,■;•, .!•, J 2 [3,3,3] 4 [4,3,2] . . 

( l 2 3,a " 1575 6 Z wsj (4, 3) + 2025 6 Z^j(4,3)' 1 j 



which simplifies to: 



(T 4 T 3 r2) = (59 + 14a + a 2 ) (4 + a) 2 a (l + a)(2 + a)(3 + a) 
2 (7 + a) 2 (8 + a) (9 + a) 2 (10 + a) (11 + a) (12 + a) 

See fig. [D /or otter examples of evaluation of (T* TfTf ^ These evaluations have taken a few seconds on a 
standard laptop. 

Example 5 Consider the same average (T\ T2T3 ) a for (3=1. The expansion of the monomial function m^^] i n 
the Jack basis is (for an alphabet of size N = 3) 

_ 1 t(2) , 1 j(2) 

TO [4,3,2] - 504[)0 J [3 , 3)3 ] + lgl44 ^[4,3,2] 



20 30 

a 



FIG. 2: Values of (T* T 2 3 T 3 2 ) a as a function of a £ [0, 50] for (3 = 4. 



After substitutions, one obtains: 



which simplifies to 



{T?TiTi) a = 



L 



(2) 

3,3,3] 



r(2) 
[4,3,2] 



(Ti 4 T 2 3 T|) 



50400 6 ^=^1,3) 18144 6 ^=7(1,3)' 

(17 + 8a + a 2 ) a (1 + a) 2 (1 + 2 a) (3 + 2 a) 
(6 + a) (5 + a) (3 + a) (4 + a) 2 (2 a + 7) (9 + 2 a) 



(39) 



(40) 




FIG. 3: Values of (T{ T 2 3 T 3 2 ) a as a function of a £ [0, 50] for (3 = 1. 



See -Fig. [H/or other examples of evaluation of (TfT^T^a. 



Asymptotic analysis for TV — > oo: The computations (see e.g. Fig. [21 and EJt a nd eq. (|D18|) suggest asymptotic 
behaviors similar to the case (3 = 2. We will detail this in a forthcoming paper [20| . 
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IV. CONCLUSIONS 



In this paper we considered nonlinear statistics of permanents on the eigenvalues of classical invariant random 
matrix ensembles. Motivated by applications to the problem of quantum transport in chaotic cavities, we first gave 
general formulae (based on a hyperdeterminant version of the Heine identity) for averages of permanents, valid for 
P = 2 , any weight function ui(x) and any set of permanent functions {ipi(x)}. The question of numerical efficiency is 
then addressed, and much quicker algorithms arc found for the specific case of Jacobi weight with 7 = 1: the analysis 
is based on symmetric functions expansions, whose main merit being that the complexity does not grow at all with N 
(the number of integration variables). This results first in a remarkable increase in efficiency, and secondly it assists 
performing an asymptotic analysis for large N. This reveals an interesting combinatorial structure lurkin g b ehind, 
and a detailed analysis of the factorization conjecture we put forward is deferred to a separate publication [201 ]. 
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APPENDIX A: QUANTUM TRANSPORT IN CHAOTIC CAVITIES: THE RANDOM SCATTERING 

MATRIX APPROACH 

Consider an open chaotic cavity of sub-micron dimensions with Ni and N 2 electronic channels in the two attached 
leads. Once the system is brought out of equilibrium by an applied voltage, it is well established that the electrical 
current flowing through such a cavity displays time-dependent fluctuations, associated with the granularity of the 
electron charge e, which persist down to zero temperature [26]. 

We consider here the Landauer-Biittiker scattering approach 0, [U HH. This amounts to relating the wave 
function coefficients of the incoming and outgoing electrons through the unitary scattering matrix 5* (2N x 2N , if 
N = N 1 +N 2 ): 

s = (; (ad 

where (t, t') and (r, r') stand for transmission and reflection submatrices among different channels. 

The theory predicts that many interesting experimental quantities are represented by linear statistics (see ([9])) on 
the eigenvalues of the N x N hermitian matrix W (if N = min(A r 1 , N 2 )): for example, the dimensionless conductance 
and the shot noise are given respectively by 5 = Tr(ttt) ^ a nd ? = Tr[tt^(l- tt^)]]MW 



m Sllil i 



Random Matrix Theory, along with insightful semiclassical approaches [3l|, [32j, [33j, [34], [35| is known to be very 
effective in describing universal fluctuation statistics in open cavities. The simplest assumption is that the scattering 
matrix S for the case of chaotic dynamics is drawn from a suitable ensemble of random unitary matrices [36l.l37ll38l.l39l]. 
Assuming then ballistic point contacts (26[, a maximum entropy approach leads the probability distribution of S to 
be uniform within the unitary group, i.e. S belongs to one of Dyson's Circular Ensembles pi l40l]. 

The unitarity constraint induces a certain joint probability density on the transmission eigenvalues {Ti} of the 
matrix tt\ from which the statistics of interesting experimental quantities could be in principle derived. This jpd is 
exactly of the Jacobi form with 7 = 1 considered throughout this paper [26|, [H, |4l[ : 

1 N 
P(Ti,..., T N ) = J] |T, - T k f J] T?- 1 (A2) 

where the Dyson index f3 characterizes different symmetry classes (/? — 1,2 according to the presence or absence of 
time-reversal symmetry and (3 — 4 in case of spin- flip symmetry), and 

a=|(|JVi-JV 3 | + l) (A3) 

The eigenvalues Tj are thus correlated random variables between and 1 and have an intuitive interpretation in 
terms of the probability that an electron gets transmitted through the ith channel. From (|A3|) . assuming N\ = £N 2 
one has a(N) - ((3/2)(i - l)N for large N (see eq. (|D3l ). 
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From (|A2|) . in principle the statistics of all the interesting quantities can be tackled, such as: 



3 = Tj (conductance) (A4) 

i=l 
N 

•J> = ^2T i (l-T i ) (shot noise) (A5) 

i=l 
N 

T p = Jf (integer moments) (A6) 



along with any other linear statistics A = Y^,iLi /(^)- 

Linear and nonlinear statistics. The average and variance of the above quantities are known, both for 
large N [TJ, [2(| [42j and, very recently, also for a fixed and finite number of channels Ni , N% [HI, EH, E3] . The 
full distribution of the quantities above has started recently to be the subject of thorough investigations: for the 



conductance, a formula was derived for N-\ = N2 = 1, 2 I45l |46|. 14711 . Until very recently, the distribution of the shot 
noise was known only for Ni — N2 = 1 |48|. Then, in 

mm 

formulas for the distribution of conductance and 
shot noise, valid at arbitrary number of open channels and for any /3, are derived, among other many interesting 
results. In and [5l|, recursion formulas for the efficient computation of conductance and shot noise cumulants 
were reported. In a recent letter Q, a large deviation approach to the problem of finding full distribution of linear 
statistics in chaotic cavities (valid for a large number of open channels in the two leads) was put forward. At odds 
with the linear statistics, about which virtually everything is known, the nonlinear statistics mainly considered in 
this paper (T^ 1 ■ ■ ■ T N N ) a is more tricky. It appears naturally when considering moments of linear statistics such as 
(S n )a or {'J' n ) a , as well as covariances of linear statistics such as cov(S, T), after expanding the above mentioned 
averages using the multinomial theorem. Results about these objects have recently appeared [HI, EH, E3 and the 
present work is yet another step in the same direction. 

APPENDIX B: HYPERDETERMINANTS 

The notion of hyperdeterminant was first defined by A. Cayley in 1843 during a lecture at the Cambridge Philo- 
sophical Society, about the possibility of extending the notion of determinant to higher dimensional arrays. The 
simplest generalization is given for a kth order tensor on an n-dimensional space M — (Mi lt .„ t i h )i<i lt .„ t i k < n as: 

Det M = — V e(a)M a (Bl) 

Cr=(CTl,...,CTfe)eSj 

where e(cr) is the product of signatures of the k permutations, M a = M ai n v... )0> (i) ■ ■ • Ma 1 (n),...,a k (n) an d &n is the 
symmetric group. It is straightforward to see that Det M = if k is odd. 

A further refinement is due to L. Gegenbauer (see e.g. [52|) in 1890 who generalized (|B1[) to the case where some 
of the indices are non-alternated. More precisely, if I denotes a subset of {1, . . . , k} one has: 

Det z M=i ]T e m ffi F ( B2 ) 
' o-=(a-i,...,CT fc )ee£ \iei / 

In particular, in the main text we defined: 

Det + (M ill<a , is ) = Det {2i3} (M llil2il3 ) (B3) 

No matter how many indices are non-alternated, every hyperdeterminant admits an expansion in sums of lower- 
order hyperdeterminants. More precisely, each hyperdeterminant of dimension k (where the dimension is just the 
number of indices) and order n is equal to a linear combination of (nl) k ~ e hyperdeterminants of dimension I and order 
p. The hyperdeterminants in the sum are obtained by fixing some of the indices. In particular, it is always possible 
to expand a Gegenbauer hyperdeterminant (see e.g. [521 ]) as a sum of (n!)*" 2 conventional determinants: 

Det 7 M= e det(Af CT3 — fffc ) (B4) 

\i£l ) 
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where M (73 '---' (7k denotes the n X n matrix such that (M (T3 >—' (Tk )i > j = M ii j :Cr3 u\..^ ak t i y 

Combining (IB4[) with (|B3|) , one easily obtains the expansion in (fTB|) . Note that a more general version of (|B4|) can 
be found in 15311 - 



APPENDIX C: PARTITIONS AND SYMMETRIC FUNCTIONS 



A partition is a finite sequence A = (Ai, A2, . . . , Ajy) of non-negative integers (called parts) such that Ai > A2 > 
. . . > An > 0. We define the weight of a partition |A| as the sum of its parts, and its length, £(X) as the number 
of its non-zero parts. Two partitions differing only by the number of their zero parts coincide. One can think of 
unidentical partitions of weight N as different ways to write the integer N as sums of positive integers. For example, 
one has only one partition A = (1) in the case of N = 1, but two partitions A = (2,0), (1, 1) for N — 2 and three 
A = (3, 0, 0), (2, 1, 0), (1, 1, 1) for TV = 3. 

The symmetric functions are polynomials in several variables X = {x±, . . . , x n , . . . } which are invariant by permu- 
tation of the variables. The set of all these polynomials for a given alphabet is an algebra A. In the case where there 
is no relation between the variables (this implies in particular that the alphabet is infinite) , the elements of the bases 
of the space A are indexed by partitions. This is the case, for instance, for the monomial functions which are defined 
by 

m A (X) = i £ (CI) 

where the already defined symbol A ! = n»ii' if A = [Ai > A2 > • • • > A&] = [. . . P* . . . 2 J ' a l 3 ' 1 ] and Xk > 0. Now, 
if we orthogonalize this basis w.r.t. the standard scalar product over the symmetric functions, we obtain the basis 
of Schur functions s\. The Gram-Schmidt algorithm allows to write monomial functions as a linear combination of 
Schur functions. For a given partition A, the Schur polynomial is defined as: 

det (x^ +n - J ) 

s\(xx,...,x n ) = t x = (C2) 

det (xr j ) 

^ ' l<i,3<n 

The denominator in (|C2[) is the Vandermonde determinant JX^^ — x j)- For partitions composed by just one 



part, A = (r), Schur functions are just the complete symmetric functions, sr r \(x) = h r 24J| . while for partitions of the 
form A = (1, . . . , 1) = (l r ), the Schur functions snr\ are the elementary symmetric functions e r (x). Schur functions 
corresponding to partitions of N form a basis in the space of homogeneous symmetric polynomials of degree N, so 
that any homogeneous symmetric polynomial can be written as a linear combination of Schur functions. 
There is an efficient way to compute such expansions. Let us first see an example: 

Example 6 Suppose that we want to compute the Schur expansion of m^ iy We consider the alphabet X = 
{x\, X2, £3, X4}. Evaluated on X, the monomial function gives 

m [3 ,i ] (X) = x 3100 +x 3010 +x 3001 +x 03W +x 0301 +x 0031 + 

2,1300 _|_ ^,1030 _|_ ^,1003 _|_ 3,0130 _|_ 3,0103 _j_ 3,0013 



'2 



C3 3 x % £ . Hence, 



TO [3,1](X) = S3100 + S3010 + S3001 + s 0310 + S0301 + S0031 + 
S1300 + S1030 + S1003 + S0130 + s 0103 + S0013 

= S3100 + + - s 2 ii + + smi+ 

-S2200 +O + S11U +0 + + 
— S31 — S211 — S22 + 2sim. 

The standard algorithm thus goes as follows (see exercise 11 pag. 110 in [241]): 
1. Expand the monomial symmetric function m\ in the variables x 1 

mx = J2 X * ( C3 ) 
where / stands for all distinct permutations of A considered as a vector of size N, completed by 0s if necessary (ex: 



[3,1] - [3,1,0,0] for N = 4), an d x 1 



. x 



Js 
N 
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2. Replace each xi by a generalized Schur function si, denned as: 

det(£/ +n ~ J )i<.j ,■<„ 

■s/ = — ^ / \ = de H s [/ ( -i+j]h<»,.j<*(i) (L4) 

lli<jV x i *jv 

where so = 1 and S-i — for each i > 0. Note that such generalized Schur function is equal to a traditional Schur 
function times a coefficient or ±1. 

3. Replace each si by 0, ±1 times the corresponding Schur function, according to the rule for i < j: 

-s...,j-M+i,... if i < j — 1 , C5 ^ 
if i = j — 1 ' 

Another important example of basis is given by the Jack polynomials which are a one-parameter deformation of the 
Schur functions. We follow the notation of [24j. One starts from the deformation of the usual scalar product defined 
on power sums by 

(px,P»h =Z iW zx5 x>li . (C6) 
where P[Ai,...,A fe ] =P\i-- -P\u and Pn = J2 x ex x " 1 ■ Tne coefficient z x is given by: 

z x =l[a i \i a \ (C7) 



di being the number of occurrences of i in A. 

A 



The Jack basis P x is obtained orthogonalizing the monomial basis with respect to the dominance order -<. This 



means: 

1. (Pf ,pP > e = 0if A^a* 

where /i -< A means J2i=i Mk — Ei=i f° r a ^ K - 
Example 7 One has 



m \i.i,i] = p < 



Cfl 

1,1,1]' 



m, , — P W i ' l^ 1 ' 1 1,1,1 P 
TO[9 11 — -frn ii + ._/») —(J) — ir 

[1,1,1] 



[2,1] — J [2,1] "T" /p(5) p<« v - 1 [1,1,11 

\Ml,l,ll'Ml,l,lK« 



- p(?) _ 6 p(g) 
^[2,1] e+2^[l,l,l] 

_ _ p (g) (^Pl-Pp,!])! p(g) , (^Mgl^l^ p(g) 

m [3] - ^[3] + (plfl P (0 \/ [2,1] + /p«> P <£) )c Ml, 1,1] 

\ [2,i] '-^[2,i]'£ ni,i,i]' r [i,i,i]'{ 

_ p«) 3_ p(0 , 6 p(S) 

- M3] 2 5+1 M2,l] + K+2)(£+l) [1,1,1] 

Unlike the Schur functions, the Jack polynomials pjp are not orthonormal. Many normalizations are encountered in 
literature. The normalization J x which is used in our paper is the integral form of Pjp: 



jP = H(Za x (s) + l x (s) + l)p( 



A 

sGA 



where the product is over the nodes s = of the partitions A (regarded as a tableaux), a\(s) = Xi — j and 

l\(s) = X'j — i if A' denotes the conjugate partition of A. 
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APPENDIX D: ASYMPTOTIC ANALYSIS FOR N -> oo AND /3 = 2 



The main formula (|35[) lends itself to a quite interesting asymptotic analysis for N — > oo. Since the sum does not 
depend on N, one may be tempted to analyze the large N asymptotics of individual summands. Quite interestingly, 
this is not sufficient: the individual summands actually diverge when TV — > oo, whereas the full /z-sum converges as it 
should. More precisely, each individual summand factorizes into the product of 0) a convergent term depending of a 
and □) a divergent term with a polynomial asymptotic behavior but with no a-dependence. Indeed, if one replaces 
the normalization constant by its explicit value, one can cast (|35j) in the form: 



M Hi AT" ' 

\i ^ A -ll j-i + l i - i - i - i -2N + n i -i-j + a K ' 

M i<j i=lj=l 

* » ' >> . ' 

(□) (0) 

where M = l(X) < N does not increase with N. 

No matter what the dependence of a on N is, the fact that (0) converges for N — ► oo is evident being the ratio of 
polynomials in N of the same order. 

We are interested in computing the following limit: 



£ A = hm (T^---T^) a{N) 

JM —>oc 



(D2) 



where the dependence of a on TV is arbitrary. In the context of the present study (see Appendix |A| , the parameter 
a is a linear function in N 

a(N) = ^(£- 1)N+ 1, (here/3 = 2). (D3) 
but one can extend the results to the cases where a is a polynomial in N , 

a(N) = (e-l)NV+J2 b mN rn , 

with p 6 Q. Note that only the highest degree part of a(N) gives contribution in the limit N — > oo. Hence, we will 
only consider the case where 

a(N) - (£- l)N p 

with p s Q. The computation of the sought asymptotics is now straightforward using symbolic softwares, as in the 
following example. 

Example 8 Suppose one has to compute £[ 2 ] = li m iV— >oo(2i )a- First, using H35\) , one obtains after simplification of 
each summand: 

2 _ 1 (l + N)(N + a)(-l + N + a) 1 (TV - 1) (-1 + N + a) (-2 + N + a) 
{ l)a 2 (2iV + a)(2iV-l + a) 2 (2 iV - 1 + a) (2 N - 2 + a) ' { ' 

Note that each individual summand does not converge for N oo as remarked above. After simplifying the full 
expressions one obtains: 

2 _ (-l + iV + a)(-3iV + 3iVa-2a + a 2 + 3iV 2 ) 

= (2 N + a) (2 N - 1 + a) (2 N - 2 + a) ' (D5) 

This expression has the same degree in N in the numerator and denominator, no matter what the dependence of a 
on N is. So the limit exists and is given by the ratio of the highest powers of N. 

Example 9 See in Fig. [^J an example showing the convergence of (TfT^T^-sN+i- 

In all cases, numerical evidences suggest the following conjecture about L\, which will be analyzed in further detail 
in a forthcoming publication (20l |. 
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FIG. 4: Values of (T 1 4 T 2 3 T 3 2 W+i as a function of TV G [0, 10 6 ] for (3 = 2. From JD6J and fD15| for p = 1 (or equivalently 
(|DlT|) '), one has lim JV _ 00 (TTTiTgjsN+i = 1253598528/6103515625 ~ 0.20539, in full agreement with the plot. 



Conjecture 1 (Factorization of limits) Let A = (Ai, . . . , Am) be a partition. One has 

U 

This means that it is always sufficient to analyze the limit in the case of partitions with one single part. 

Example 10 Consider again the case where A = [4,3,2], but for a ~ (£—l)N. After a brief computation one obtains: 

4 3 2 (l + 3£ + 9£ 3 + 9£ 2 + 3£ 5 + 9£ 4 + £ 6 ) (l + l + i 2 ) (£ 4 + 2£ 3 + 4£ 2 + 2£+l)£ 3 
lunJ Tl T 2 T 3 ) (e „ 1)N = ^—^ . (D7) 

But one has also 

. 4 . £(l + 3£ + 9£ 3 + 9£ 2 + 3£ 5 + 9£ 4 + £ 6 ) 

lim (T-, )(e-i)N = ? 

i;jv (l + £) 7 

. J?(f 1 + 2f 5 + 4<? 2 + 2<?+l) 
i& <Tl ] ^- 1)N = (1^8 

, ^(l+^ + £ 2 ) 

lim (T-, }t£-i)N = 5 — • 

n^™ 1/(1 1)N (l + £) 3 



Hence, 



lim {T 4 T 3 T 2 ) {e _ 1)N = lim {T 4 ) (e _ 1)N lim (T?) (/ _i )JV Jim (I?) ( /_i )Ar . (D8) 



Assuming Conjecture |T]), it remains to consider the limit 

%.] = lim I [fc] (a,/V) (D9) 



N 

where 



J w (a,JV) = <2i>a (D10) 
i.e. the case when the partition A is composed by just one part A = [k]. 
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For a(N) ~ § (£ - 1)N + 1, the limit (|D9| has been computed by Novaes [Tj as: 

%1 . ,» (f (< _ = le+l) ±tip(^ (_?_,)■ (D11 , 

Guided by the numerics, we have found an equivalent expression (see (|D15[) ). whose direct combinatorial proof will 
be announced in a separate publication [20j] . In the case A = [k], the monomial function TOrM = ^2 i T^ is the power 

sum pk and the coefficient K\k],ft are well known (see e.g. 

fe-i 

m W = ^2(- 1 ) ts l(k-i),^]- (D12) 

i=0 

Plugging these coefficients in (|55|) , one recognizes, after simplification, a hypergeometric function: 

- _ T(2 N + a - l)r(— 1 + N + a + k)(N + k - I)] /2-N-a, -2N - a + 2 - k,-k + l,l- N 

[k] ^' ' ~ k\NW{-l + N + a)T(2 N + a - 1 + k) 4 3 \-k + l-N,-2N + 2-a,-N-a + 2-k' 

(D13) 

where 

/oi, . . . ,a p \ _ ( ai )i . . . {a p )j x i 

if (x)i = x(x + 1) . . . (x + i — 1) denotes the rising factorial. 

Suppose now that a ~ (£ — 1)7V P . In this case, numerical evidences suggest the following alternative representation 
for (fDTTT) : 

(wif^E,^" (Si) (f:])' 

£ w = Jim /„((<- l)iV»lV) = (^1 f0Ip<1 (D15) 

[ 1 for p > 1 

where [w] (resp. kvj) denotes the smallest (resp. largest) integer larger (resp. smaller) or equal to lo. 
Note that: 

• The coefficient ^ i *j ^ ( J has a very interesting combinatorial interpretation, since it is also the number of 

symmetrical Dyck pat hs with even semi- length 2k and exactly i peaks [54|. We will explore this property in a 
forthcoming paper [20| , where a formal proof of the equivalence between (ID15I) and (|D1 1|) based on the creative 
telescoping method will be provided. 

• Eq. (|D15|) for p = 1 can be equally well restated in terms of hypergeometric functions as: 

lirn^ t [k] ((£ - 1)N, N) = ( ^ + ^ 2fc _ 1 ( 2 F 1 (-k, -k; + £k 2 F 1 (1 - k, -k; 2; £ 2 )) . (D16) 

• The second and third case in (|D15|) are obtained from the p = 1 case upon setting £ = 1 and £ — > oo respectively. 

• From the factorization conjecture and the independence of Novaes' limit (|D1 1|) on the exponent j3 of the Vander- 
monde, one obtains the following result. Suppose (3 > and set 

N 



One has: 



-i— - - f m ■ ■ ■ dT N ■ ■ -T^ J] \Tj -T k ff[ Tr 1 . (D17) 
MN) J m » " f = \ 



lim h[%(l- 1)N, N;(3)= L\. (D18) 



JV^oo V 2 
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